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Abstract. We present results of our pulsar population synthesis of normal pulsars from the Galactic 
disk using our previously developed computer code. On the same footing, we use slot gap and 
outer gap models for gamma-ray emission from normal pulsars to obtain statistics of radio-loud 
and radio-quiet gamma-ray pulsars. From recently improved understanding of HII and star forming 
regions in the Galaxy, we develop a new surface density model of the birth location of neutron 
stars. We explore models of neutron star evolution with magnetic field-decay, and with different 
initial period and magnetic field distributions. We present preliminary results including simulated 
population statistics that are compared with recent detections by Fermi of normal, isolated pulsars. 
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INTRODUCTION 

NASA's new flagship, the Fermi Gamma-Ray Space Telescope (Fermi) launched in June 
2008, has opened a new era in y-ray astronomy, improving dramatically our understand- 
ing of y-ray emission from pulsars. Fermi has discovered over 70 normal, isolated and 
millisecond pulsars, many through blind period searches, superseding the y-ray pulsar 
database of six provided by its predecessor, the EGRET instrument aboard the Comp- 
ton Gamma-Ray Observatory. A group of 40 normal, isolated pulsars that are publicly 
available are used in this study. 

BIRTH DISTRIBUTIONS 

Previous to this study, we have used [2] the exponential, spatial distributions from the 
work of Paczyhski [5]. Our preliminary statistics of y-ray pulsars indicated that our 
simulations required a larger group of nearby pulsars, as the code was not simulating 
sufficient Fermi pulsars. As a result, we developed a new surface density distribution 
of neutron star birth locations adopting the recent discoveries of nearby H II regions 
[1]. We display in Figure 1 a comparison of the smooth exponential probability radial 
distribution [5] represented by a dashed curve, and our new model shown as a solid 
curve. This new radial distribution together with the exponential height distribution and 
a uniform angular distribution are used to generate randomly distributed neutron stars in 
the Galaxy at birth. We implement the supernova kick velocity distribution by Hobbs et 
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FIGURE 1. Radial birth probability as a function of the Galactic radius in kpc (cylindrical coordinates). 
The dashed curve represents the probability distribution given by the study of Paczyhski [5], while the 
solid curve represents our new model probability distribution implementing the newly discovered HII 
regions in the study of Bania et al. [1]. 



al. [3], and evolve the neutron stars in the Galactic potential [5] using a a Cash-Karp fifth 
order Runge-Kutta routine [6]. The birth rate of 2 neutron stars per century [4] back to 
1 Gyr in the past is assumed to be constant. 

We assume a Gaussian distribution for the initial period distribution of the neutron 
stars with a mean P of 50 ms and a width o of 50 ms in the form: 

P(P ) - e -P-Po)W (1) 

The initial magnetic field distribution is represented by a sum of two log-normal distri- 
butions as 

P(B) - ^ A . e -(rt-log(s)) 2 /^ (2 ) 

i=\ 

with the parameters given in Table 1 . In order to improve the agreement between the 
simulated pulsars and those detected by Fermi in the P — P diagram, we had to change the 
widths and means the initial magnetic field and period distributions from our previous 
work in order to produce a greater number of younger, high E pulsars. We assume 
a field decay constant of 2.8 Myr that is similar to our previous studies. We are not 
necessarily implying the existence of field decay, rather we suggest that field decay may 
be mimicking some form of pulsar spin-down. 

SIMULATION AND RESULTS 

In the process of improving the agreement with the Fermi detected pulsars, we find 
that the comparisons with the radio pulsar distributions have had to be compromised to 
some extent as can be seen in Figure 2. We represent the distributions of the indicated 
characteristics of the detected radio pulsars in solid histograms while the distributions 
of simulated pulsars are indicated in open histograms. The new radial birth distribution 



TABLE 1. Assumed initial parameters of previous [2] and current 
studies for the birth magnetic field (log Normal) and period (Gaussian) 
distributions. 
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Initial Period P a 

Previous Study [2] 300 ms 300 ms 

Current Study 50 ms 50 ms 




shown in Figure 1 leads to an excess of radio pulsars between the Sun and the Galactic 
center as shown in the comparisons of the distance histograms in Figure 2. We also 
simulate an excess of young radio pulsars as indicated in the period, period derivative 
and age histograms. 

In Figure 3 (left panel), we present the detected radio pulsars (solid dots) observed 
from a select group of ten radio surveys, and the detected Fermi radio-loud (bowties) 
and radio-weak (hourglasses), where the definition of radio-weak is used to indicate 
that these pulsars were not detected in our select group of radio surveys. In the right 
panel, we show the simulated radio and y-ray pulsars. The y-ray emission used here is 
from the polar cap/Slot Gap (SG) model. We normalize the simulation to the number 
of detected radio pulsars observed in the ten radio surveys. The radio luminosity is 
adjusted to produce a neutron star birth rate of 2.1 per century [4]. Further details of 
the y-ray beam geometry and luminosity as well as alternate y-ray models are presented 
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FIGURE 3. Period derivative - period plot of detected (left) and simulated (right) pulsars. Radio pulsars 
indicated in solid dots, Fermi radio-loud (bowties) and radio-weak (hourglasses) pulsars. The dashed 
lines represent constant magnetic field loci assuming vacuum, dipole spin-down, and dotted lines indicate 
curvature radiation and non-resonant Compton Scattering pair death lines. 



in the companion paper of this conference proceedings by Pierbattista et al. Here we 
are using a group of 30 radio- weak and 10 radio-loud Fermi pulsars that have been 
made public. Our simulation predicts a slightly smaller ratio of radio-weak to radio-loud 
with 44 radio-loud and 17 radio-weak y-ray pulsars. Preliminary studies suggest that 
it should be possible to simulate both Fermi and radio pulsars with the same group of 
evolved neutron stars. However, further refinements are required in order to improve the 
agreement between detected and simulated distributions. 
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